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BUCKLING OF SEGMENTS OF TOROIDAL SHELLS 


By Manuel Stein and John A. McElman 
NASA Langley Research Center 


SUMMARY 



> 


Nonlinear differential equations of equilibrium and buckling equations are derived 
for segments of toroidal shells near the equator and for segments near the crown. The 
equations are derived for shallow shell segments by including appropriate prescribed 
initial displacements in the nonline ar, flat plate, strain-displacement equations and by 
varying the total potential energy of the system. Closed form solutions to the buckling 
equations are obtained for simply supported segments near the equator having either posi- 
tive or negative Gaussian curvature under pressure loading with various inplane support 
conditions. Results are presented in the form of charts showing buckling coefficients as 
a function of a curvature parameter associated with the girth of the shell and a param- 
eter associated with the ratio of principal curvatures. In many instances the results 
indicate significant deviations in buckling stress for the toroidal shells over the 
buckling stress for the corresponding circular cylindrical shell under similar loading 
and support conditions. 

INTRODUCTION 



Shells of double curvature are common in aerospace vehicle structures, and buckling 
is often an important design consideration for such shells. In this paper nonlinear dif- 
ferential equations of equilibrium and buckling equations are derived for segments of 
toroidal shells, a type of double curvature shell which is attracting considerable 
interest at the present time. Solutions to the buckling equations are obtained for seg- 
ments of toroidal shells near the equator having either positive or negative Gaussian 
curvature (see fig. l) subjected to various pressure loadings. 

The nonlinear equilibrium equations are derived for shallow shell segments by 
including appropriate prescribed initial displacements in the nonlinear, flat plate, 

strain-displacement relations and by varying the total potential energy of the system. 
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These equations reduce to the large deflection Donnell equations for the case of a cir- 
cular cylindrical shell, and to the Marguerre large deflection equations for the case of 
a shallow spherical cap. The buckling equations are derived from the nonlinear equations 
in a rigorous manner by assuming the changes which occur at buckling to be small. One 
set of equations is obtained which is applicable to segments of a toroidal shell near the 
crown, and another set is obtained which is applicable to segments near the equator. For 
segments near the equator, the classical assumption of constant deflection prior to 
buckling leads to buckling equations which are the same as those given by Becker in ref- 
erence 1 for shells of double curvature having constant, but not necessarily equal, 
principal curvatures. 

Closed form solutions are presented to the buckling equations (obtained by using the 
classical assumption) for segments of toroidal shells near the equator subjected to 
lateral pressure. The assumed edge support conditions are simple support either with 
zero edge displacement, with hydrostatic pressure loaded edges or with freedom for over- 
all edge extension in the axial direction. Results are presented in the form of charts 
showing buckling stress coefficient as a function of a curvature parameter associated 
with the girth of the shell and a parameter associated with the ratio of principal curva- 
tures. The results indicate significant deviations in buckling stress for the toroidal 
shells over the buckling stress for the corresponding circular cylindrical shell under 
similar loading and support conditions. 

For many buckling problems involving deep shells, shallow shell analysis should give 
engineering estimates to overall buckling loads. In the present paper an estimate of the 
external buckling pressure for a complete torus is obtained on the basis of the study of 
the shallow segment near the outer equator. 


a 

A,B, C 
D 


SYMBOLS 

radius of curvature (figs. 1 and 2) 
constants 


flexural stiffness of shell w all . 


Et 3 

12(1 - H 2 ) 

2 - 




% 1 



E 


Young’s modulus 


k 

l 


buckling coefficient, 
length of shell 


pr l 


2 


Dn 2 


m, n integers 

NxjNyjNxy stress resultants in rectangular coordinates 


stress resultants in cylindrical coordinates 

p lateral pressure - positive in positive w direction 

p,0,z cylindrical coordinates 


r 

R 

t 


u,v,v 

x,y,z 


Z 


radius of shell equator (fig. l) 

central radius of toroidal segment near crown (fig. 2) 
shell wall thickness 

displacements, tangential and normal to the shell neutral surface 
displacements in x and y directions 
initial deflection 
rectangular coordinates 


curvature parameter. 




a 

3 

€ x'V 7 *y 
€ p> € e^pe 


curvature parameter, r/a 

buckle wavelength parameter, nZ/itr 

direct strains and shearing strain in rectangular coordinates 
direct strains and shearing strain in cylindrical coordinates 


[i Poisson’s ratio 

H total potential energy of shell 


= V 2 ? 2 


where 


V 2 


is the Laplacian operator in two dimensions 
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Subscripts: 

A refers to prebuckling displacements and stress resultants 

B refers to buckling displacements and stress resultants 

A comma indicates partial differentiation with respect to the subscripts following 
the comma. 


NONLINEAR DIFFERENTIAL EQUATIONS OF EQUILIBRIUM 

In this section the nonlinear differential equations of equilibrium are derived for 
shallow segments of a torus near the equators and near the crown. For segments near the 
equators the equations are derived in the rectangular coordinates of an osculating plane. 
For segments near the crown the equations are derived in plane polar coordinates. In 
both cases the equations are derived from the strains of nonlinear flat plate theory 
including initial deflections (see ref. 2) using the minimum potential energy method to 
obtain equations of equilibrium by the application of a variational procedure. 


Segments Near the Equators 

The strains for thin flat plates with u and v the displacements in the x and 
y directions, respectively, and with initial deflection w 0 and additional deflection 
w are given in rectangular coordinates as follows: 


€ x = V + |( w ,x) 2 + «,x w o,x - ZW ,XX 

£ y = v ,y + §( w ,y ) 2 + w >y w °>y " zw ^yy 

7xy = »y + \x + » x»,y + *,x*o,y + *,y*o,x - 




2zw ,xy) 


( 1 ) 


For an initial deflection corresponding to segments at the outer equator of a torus of 
radius r, which has positive Gaussian curvature and which is taken here to have constant 
curvature 1/a in the meridional direction (see fig. 1), w Q is taken to be 


w o 


y£ _ x£ 
2r " 2a 


(2a) 


% 
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The assumption that the initial deflections can be represented in quadratic form is con- 
sistent with shallow shell approximations. For segments at the inner equator the initial 
deflection has negative Gaussian curvature and the corresponding equation for w 0 is: 


*o 



/ r>n.\ 

K ^ ) 


If equation (2a) or (2b) is substituted in equations (l) and if the following defi- 


nitions are used 


u 


v 



(3) 


then 


e x « u X ± I + |(v ) 2 - XX 

e y = v >y + r + |( w ,y f ~ zw ,yy } 


00 


7xy = u ,y + v ,x + w ,x v ,y “ 2zw ,xyJ 

where the new u and v can be identified as the tangential displacements of the neu- 
tral surface of the shallow shell and w can now be regarded as the normal displacement. 
In equations (3) and (4) and in the equations that follow, the convention is used that 
when there is a double sign the upper sign applies to the shell with positive Gaussian 
curvature and the lower to the shell with negative Gaussian curvature. 

The total potential energy of such shells subject to a lateral pressure p is 


11 = i fT -» iZA 6 * 2 + £y2 + 2 ^ x£y + Hr* - p ff w & (5) 

Integration in the z direction and variation according to the minimum potential energy- 
method by the calculus of variations leads to the nonlinear differential equations of 
equilibrium and consistent boundary conditions. The equations of equilibrium are: 
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where 


N x,x + N xy,y = 0 
%,y + N xy,x = 0 

V TiT AT. 

± ~ + r " ( N x w ,xx + % w ,yy + 2N xy w ,xy^ = P 


(6) 



( 7 ) 


J 

Note that these equations with a -» oo are the Donnell large deflection equations for a 
cylinder and, if in addition r -> oo^ then they reduce , of course, to the von Karman large 
deflection equations for a flat plate. 


Segments Near the Crown 

The strains for a thin flat plate with an initial deflection w 0 are given in polar 
coordinates as follows: 


e p = u ,p + |( w ,p) 2 + w , p w o, P - 


PP 


= iy, + - + -i-(w n ) 2 + -if- W nW n A-zf-^-Wan+i-W 

0 P >9 p 2p 2 ' > 6 ' P 2 \p 2 P * { 

7 n = i u n + V n -— +iw_Wn+iw 0 W n Q + i w fi W- _ - 2Z W - “T Wo 
'p0 P ,0 jP P p ;P P jP P °>P ip >P y p2 


(8) 


where w is the additional deflection. For an initial deflection corresponding to seg- 
ments at the crown of a torus of central radius R where the torus is taken here to have 
constant meridional curvature l/a (see fig. 2), w 0 is taken to be 


w 0 = . (p - ._ R d (9) 

2a 

As before, the assumption that the initial deflections can be represented in quadratic 

form is consistent with the shallow shell approximation. Thus 
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( 10 ) 



The total potential energy of such a shell subject to a lateral pressure p is 

= v,) fff{ £ p 2 + € 0 2 + 2 ^ e p € e + 7 pe) p d P de dz - P JJ W P d P de ( dl ) 


Integration in the z direction and variation according to the minimum potential energy 
method by the calculus of variations leads to the nonlinear differential equations of 
equilibrium and consistent boundary conditions. The equations of equilibrium are 


Np,p + i(N p - N 0 + N pe?0 ) = 0 

|( N e,e + ^pe) + N p e, p = 0 


/ (12) 


ikAw + ^ + £ " R 


a a N e - (V,PP + \ N 0 W ,00 + V,p + I N pe w ,p0 - \ N p0 w ,0j = p 

v p K p / j 


where 


N P = 


Et 


( € P + 


z=0 


1 - H 




N 


pG 


2(1 + p)( 7 P 0 )|z=C 


(13) 


These equations with R = 0 reduce to the Marguerre large deflection shallow spherical 
cap equations. 


BUCKLING EQUATIONS 

Buckling equations are derived on the premise of bifurcation behavior with the non- 
linear equations just derived used to determine the deformations and stresses prior to 
buckling a.nd to determine the buckling equations. 
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For the problems considered, the shell prior to buckling deforms axi symmetrically . 
Small changes from this configuration, not necessarily axisymmetric changes, are con- 
sidered in order to obtain the buckling equations. Buckling equations are also derived 
using the assumption that the deflection w is constant prior to buckling. 


Segments Near the Equator 

Prior to buckling for the problems considered the deformations would be axisym- 
metric; thus equations (6), with the variables functions of x only, would apply. Thus, 
for axisymmetric deformations 


where 


N xA,x = 0 
NxyA,x = 0 

± ^ + 


Dw 


A, xxxx 


n xA w A, 


XX 


= p 


(14) 


n xA ~ 


Et 


1 - 


I wa -i o 

U A, x 1 T + 2 W A, x 


NyA = 


Et 

4k 

1 - n 2 

r 


4 * -4*. * 4 


Et 


NxyA = 2(1 + n) V A,x 


+ 




*A 

r 


1 

2 



To the prebuckling deformations obtainable from these equations and the boundary condi- 
tions, small changes that occur during buckling may be added ^u = u A + u B , v = v A + v B , 
w and the sum should also satisfy equations (6). The above relations may be 

subtracted out after the sums are inserted in equations (6), and terms of higher degree 
than linear in the buckling displacement may be neglected to give the buckling equations 
which follow: 
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( 15 ) 


h xB,x + N xyB,y “ 0 
%B,y + N xyB, x = 0 

■°^ W B ~ a + ( N xA w B,xx + N xB w A, xx + ^yA'fyyy + N yB w A,yy 

+ 2NxyA w B, xy + ^xyB^xy) = 0 

where 


NxB = 


Et 


MyB = 


i - n 2 L 


Et 


%* 1 T + Xx W B,x + ^Ky + T + W A,y W B,y 


1 - n‘ 

Et 


V B,y + T + W A,y W B,y + % U B,x ± — + W A,X W B,X 
?xyB = 2 '(l + n) ( U ^y + V B,x + w A,x w B,y + w A,y w B,x) 


These equations have variable coefficients and would be quite difficult to solve for many 
cases . 

If, instead, the assumption analogous to that of classical cylinder buckling theory 
is made, that the deflection w is constant prior to buckling (w A = Constant then with 
p, the pressure, a known constant it may be seen from equations (lV) that and N^ A 

must be constant also. Equations from which these constants may be determined for differ- 
ent inplane end conditions are 


w Et / , V A . w a\ 

NxA ~ n 2 \ A, x “a ^ r / 
1 - \± 

•Va - - 4a,x * ^)] 




(16) 


+ !!** + 5[A = 

a r F - 

In these expressions u A ^ x is also constant, and therefore, for the special case of 

u A = 0 at both ends u A ^ x would be zero. Three different inplane end conditions are 

considered for the problems solved in this paper and are discussed in a subsequent 

section. 
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With w^ a constant the buckling equations can be written 


where 


N xB,x + N xyS,y " 0 

%B,y + N xyB,x = 0 l 

ri ^ w B 1 "JT + -%r- - ( N xA w B,xx + N yA w B,yy + ^xyA^xy) = ° 


w xB 


Et 


Wp 

Up v + -^ + 

b, x a 


N yB = 





^B.y 

n( u B,x 



Et / \ 

' IuTToK* + v *,*) 


(17) 


These buckling equations, obtained through the classical assumption, have constant coef- 
ficients and agree with those derived by Becker in reference 1; with a -> oo they agree 
with the Donnell equations for buckling of a cylinder. 


Segments Near the Crown 

Prior to buckling, for the problems considered, the deformations would be axisym- 
metric; thus equations (12) with the variables functions of p, only, would apply: 


Np,Ap + £(N P A - N 0A ) = 0 


§ N p6A + N p0A,p = 0 


D 

P 





a 



(' 


pA w A,pp 



= P 


where 


N pA 
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N0A 


Et 


1 - 


IX/ , 1 

2 

p - R \ 

T + TA,P + 2 

W A 

A,p 

a Xp) 


P® A 2(1 + (i) 


_Et / _ VA\ 

L + u) IX P P / 


Proceeding as for segments near the equator, the buckling equations are 

N pB,p + p( N pB “ N 0B + N p9B,e) = 0 

^( N 0B,e + ^pgb) + N P0B,P = 0 

+ “S“ + N 0B " ^ N pA%PP + N pB w A,pp + ~2 N 9A w B,ee + ^2 N 0B w A,00 + N 0A w B,p 

+ jj N 0B w A,p + f N p0A w B,p0 + f N p0B w A,p0 “ ^ N p0A W B,0 " ^ N p0B w A,0 ' = 0 

where 


NpB = 
n 0B = 
N, 


Et 


1 - P 2 
Et 


..21o 

- H 


Et 


u B,p + w A,p%p - ( £ -T^) W B,p + ^ V B,0 + T + ^ W A,0 W B,0j 

,p w b,p - (^it) w b,p 


(l U B 1 

fe Xe + - + ^ XeXe + » 


Xp + X 


1 

V B 

1 

, 1 

(p ~ I 


p U B,0 

+ V-R „ 

B >P p 

+ p W A,p w B,0 

+ p w B,p w A,0 " 

V pa 

7 B,0 


P® B 2(1 + n) 

For the segment of the toroidal shell near the crown, it is expected that the classical 
assumption of constant deflection w^ prior to buckling will give reasonable stress 
resultants. With this assumption the prebuckling stress resultants can be found from 


NpA,p + i(N pA - N 0 a) = 


p N p6A + N p9A,p “ 0 
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N 0 A = P 


N pA + p - R 
a pa 


The buckling equations then become 


where 


Vb,P + |( N pB - N 6 B 

K N eB,e + ^pgb) 


+ £_i_ 


pa 


■ n 6 B 


+ | N 9 A w B,p + § N p 0 A w B,p 0 


+ N p 9 B,e) = 0 
+ N p 9 B,p = 0 
N pA w B,pp + ^2 N 9 A W B ,99 
jjj N p 0 A w B, 9 ^ = 0 


pB = iT7p B ' p ■ (“H^Kp + 4 t b,( 


Et a u B 

n ob = — zf v B,e + - + P 


: u b,p - (^) w 


Np0 B = — 


Et 


(1 + n) 


1 V B / 

p u B,e + V B,P - — ~ \ 


V B (p - R 


pa 
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SOLUTIONS FOR SEGMENTS NEAR EQUATORS 

Closed form solutions to the buckling equations just derived are now presented for 
simply supported segments of toroidal shells near the equators under pressure loading 
with various inplane support conditions. The shells considered extend completely around 
the equator, and the equator lies at their midlength (see fig. l). Shells of both posi- 
tive and negative Gaussian curvature are considered. 

The inplane support conditions considered are listed below together with the N^ 
and NyA determined (on the basis of the classical assumption that w a is a constant) 
from equation (l6): 
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(a) zero edge load 


^xA - 0 
N yA = pr 


(b) hydrostatic pressure loaded 



With no applied shear stress, ~ °* 

The buckling equations ( 17 ), obtained through the classical assumption may be 
written in terms of u-g, v s , w-g as (dropping the subscript B) 



With the origin now taken along the lower edge, the simple support boundary condi- 
tions on the buckling displacements at x = 0 , l are 

w=w xx = v=N x = ° (20) 
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The buckling pressure coefficients have been calculated for Poisson’s ratio 
M- = l/3 anc ^ the results of the calculations are plotted in figures 3 to 6. For the 
shell of positive Gaussian curvature under external pressure the results for higher 
values of Z lie along a straight line as shown on the logarithmic plot. Simple 
results for the critical pressure obtained analytically (assuming 3 large) for any 
value of Poisson’s ratio are presented below for these straight line regions: 


- Ik - 


-p = 


[5(1 - n 2 ) 


Et£ 

1/2 ar 


for zero edge loading 


-P = 


-p = 




Et g 

ar 


for hydrostatic pressure loading 


L r 1 'j 

1 + 2|ia + a, 2 


Et" 


(1 + a) [3(1 - n 2 )] 


# for zero ed«e displacement 

1/2 ar “ 


ESTIMATE OF EXTERNAL BUCKLING PRESSURE FOR COMPLETE TORUS 


Prebuckling values for the stress resultants for the pressure loading of a complete 
torus analogous to those obtained for the shallow shell (see eq. (20)) are 


n xA 


p»(i - 5)'' 


%-f J 

and proceeding with these values as for the shallow shell leads to 


(23) 


k = 


(1 ♦ ♦ if (1 ♦ 


(24) 


Finally, for p large, the buckling pressure can be estimated to be (probably not 
valid for a circular torus when r/a is much greater than 2) 

k2 


-P * 


Et £ 


[3(1 - 4 : 1/2 * 2 


(25) 


This result indicates that the torus buckles when the circumferential stress reaches 

0.6 — . The classical results show that a cylinder of radius a in axial compression 
a 

buckles when the axial stress reaches 0.6 =— and the circumferential stress is zero. 

SI 

The sphere of radius a under hydrostatic pressure buckles when the stress (in any 

direction) reaches 0.6 — . However, the torus under hydrostatic pressure buckles when 

a 
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the circumferential stress reaches 0.6 — even though the meridional stress exceeds 

8L 

this value. 

DISCUSSION OF RESULTS 

Nonlinear equations for segments of toroidal shells near the equator and near the 
crovn have been derived from flat plate equations by including appropriate prescribed 
initial deflections. Buckling equations have been derived for both kinds of segments. 
Utilizing the classical assumption of constant deflections prior to buckling, buckling 
equations with constant coefficients have been obtained for segments near the equator. 
These equations have been solved in closed form for pressure loading of simply supported 
segments with both positive and negative Gaussian curvature having three different 
inplane edge conditions: (l) zero edge load, (2) hydrostatic-pressure-loaded, and 

(3) zero edge displacement. 

Results obtained for the case of external pressure buckling with zero edge load 
(lateral pressure) are presented In figure 3* For a given value of the curvature param- 
eter Z associated with the girth of the shell the external pressure required for 
buckling increases significantly over the buckling pressure for the cylinder (r/a = 0) 
as the curvature in the meridional direction Is increased to form a shell of positive 
Gaussian curvature. The corresponding external pressure decreases significantly as the 
curvature in the meridional direction is increased to form a shell of negative Gaussian 
curvature. For higher values of Z, the curve for the cylinder (r/a = 0) has a slope 
one-half. For the spherical segment (r/a = l) the curve has a slope of unity. And for 
the corresponding shell (r/a = l) of negative Gaussian curvature the curve has a slope 
of zero; thus the curvature contributes little to the strength. Both the shell of 
positive and the shell of negative curvature do not buckle under internal pressure for 
this inplane edge condition. 

Results obtained for the case of external hydrostatic pressure buckling are pre- 
sented in figure 4. These results follow trends similar to the lateral pressure results 
and again both the shell of positive and negative Gaussian curvature do not buckle under 
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internal pressure. Again significant increases in external pressure required for 
buckling are available for shells of positive Gaussian curvature over cylinders of the 
same curvature parameter Z , and significant decreases in buckling pressure occur for 
shells of negative Gaussian curvature over cylinders of the same Z. 

Results obtained for the cases or external ana internal pressure buckling, respec- 
tively, of shells with zero edge displacement are presented in figures 5 and 6 . In the 
case of external pressure as shown in figure 5 significant increases in the pressure 
required for buckling are available for shells of both positive and negative Gaussian 
curvature over that for a cylinder at the same value of the curvature parameter Z. 

Due to the u = 0 condition, tensile stresses develop at the edges that tend to sta- 
bilize the shell of negative curvature as the ratio of radii r/a increases. The 
shell with negative Gaussian curvature also buckles under internal pressure with this 
zero edge displacement condition (see fig. 6) provided the ratio r/a is greater than 
Poisson’s ratio fi. As r/a increases for given Z compressive stresses at the edges 
increase and cause a decrease in the buckling pressure. 

CONCLUDING REMARKS 

The present analysis starts with accepted (von Karman) nonlinear flat plate strains 
including initial deflections and in a consistent and straightforward manner derives 
nonlinear shallow shell equations from which buckling equations are determined. Non- 
linear equations and buckling equations for other shallow shells such as conical frustums 
away from the apex and segments of toroidal shells away from the equator or crown may be 
derived In a similar manner. 

Buckling pressures have been obtained in chart form for simply supported toroidal 
segments near the equators having both positive and negative Gaussian curvature. Three 
inplane edge conditions are considered. Limiting values of the buckling pressure for 
shells of large Z having positive Gaussian curvature are given in equation form. An 
estimate based on this shallow shell theory of the buckling pressure of a complete torus 
is also given in equation form. 
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The present theory Is limited by the shallow shell approximation which, however, 
may not be too serious for deeper shell buckling analysis. This conjecture that shallow 
shell theory may be used to analyze some deep shells is based on the consideration that 
shells buckle first where the curvature is most shallow. This consideration was used in 
estimating the buckling pressure for the complete torus. 

Another limitation of the present results exists because of the disagreement 
between some shell buckling solutions for perfect shells and experiment. For the cylin- 
der under external pressure, there is good agreement between theory and experiment. How- 
ever, for the spherical segment it is expected that agreement between theory and experi- 
ment will not be nearly as good as for the cylinder judging from buckling results for the 
spherical cap. Similar limitations probably arise for the buckling of other toroidal 
segments and for the complete torus. Thus, for design purposes the present results may 
only serve as a guide by specifying the buckling pressure for a perfect shell. 
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Figure 1.- Configuration for segments near equators. 
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Figure 5.- Buckling of toroidal segments with zero edge displacement under external 

lateral pressure. 





